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Abstract 

The light hadron spectrum as computed in nonperturbative QCD 
is reviewed and compared to lattice data and experiment. The mixing 
of mesons, hybrids and glueballs is calculated in the Field Correlator 
(T) . Method. The strong decay mechanisms are found out in the method 

and compared to the known phenomenological models. The role of 
sea quarks (unquenched approximation) is studied analytically using 
radially excited mesons as an example, and compared to experiment. 



1 Introduction 

> ■ 

^ ■ In the last decade there has been a substantial progress in understanding 



the QCD spectrum, both in analytical methods [T], and lattice calculations 
On the analytic side the most economic and promising turns out the 
Field Correlator Method (FCM) j^j, which starts from the basic principles 
and correct description of the QCD vacuum with the help of gauge-invariant 
field correlators. It was proved [3] that confinement can be described by 
the lowest, quadratic correlator (F IJiU (x)Fxa(y)) which contains two scalar 
formfactors D(x — y) and D\(x — y). Lattice data [H] show that D and D\ 
describe QQ static potential with a few percent accuracy, and the string 
tension a is obtained directly from D(x): a — \ J D(x)d 2 x. 

In addition D(x) = D(x/T g ) contains another important parameter - the 
gluon correlation length T g , which was found on the lattice and analyti- 
cally 7j to be very small, T g < 0.2 fm. 
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This circumstance allows to develop the local Hamiltonian and Lagrangian 
methods for the description of the qq and 3q bound states, which will be dis- 
cussed in the next section. 

However, this is not the whole story since there are valence gluons which 
can also form bound states by themselves (glueballs) and with the qq (hy- 
brids). To introduce these states one should define the valence gluons in con- 
trast to the nonperturbative background, and one can do it unambigiously 
in the framework of the background perturbation theory jHj. In this way one 
gets the local Hamiltonian also for hybrids [3 [TUJ and for glueballs [IT] and 
calculate the corresponding spectra in good agreement with lattice data. In 
doing so one realizes that a gluon excitation "costs" around 1 GeV, which 
allows to disregard these states in the first approximation when computing 
the lowest meson or baryon states. 

However the exact treatment requires the introduction of Fock tower of 
states, and consideration of mixing between meson and hybrid states, which 
is done in Section 3. The effects of sea quarks on the spectrum and strong 
decays is considered in section 4, while Section 5 is devoted to conclusions. 

2 Hamiltonian 

There are two possible approaches to incorporating nonperturbative field 
correlators in the quark- ant iquark (or 3q) dynamics. The first has to deal 
with the effective nonlocal quark Lagrangian containing field correlators ^2] • 
From this one obtains first-order Dirac-type integro-differential equations for 
heavy- light mesons J21EI]> light mesons and baryons j!4j . These equations 
contain the effect of chiral symmetry breaking [12, which is directly connected 
to confinement. 

The second approach is based on the effective Hamiltonian for any gauge- 
invariant quark-gluon system. In the limit T g — > this Hamiltoniam is simple 
and local, and in most cases when spin interaction can be considered as a 
perturbation one obtains results for the spectra in an analytic form, which 
is transparent. 

For this reason we choose below the second, Hamiltonian approach ^5] 
Hp] . We start with the exact Fock-Feynman-Schwinger Representation for 
the qq Green's function (for a review see ^7]), taking for simplicity nonzero 
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flavor case 



p OO P OO 

G% y) = j o d Sl J ds 2 (Dz) xy (Dz) xy e- K ^x 



(trT in ( mi - £> 1 )W CT (C)r 0Ut (m 2 - D 2 )) A (1) 

where K { = J Sl dr^mi + ^(i^ } ) 2 ), T in)OUt = 1,75, ... are meson vertices, and 
W a (C) is the Wilson loop with spin insertions, taken along the contour C 
formed by paths (Dz) xy and (Dz) xy , 

W a {C) = P F P A exp{tg f A^dz^x 

PS 1 C '52 

xexp{g a^F^dn-g a { ^F^dT 2 ). (2) 

•J " 

The last factor in (J2J) defines the spin interaction of quark and antiquark. 
The average {Wo) a in (0) can be computed exactly through field correlators 
(F(1)...F(ti))a, and keeping only the lowest one,(F(l)F(2)), which yields 
according to lattice calculation |S] accuracy around 1% .4 , one obtains 

(W a (C)) A ~exp(-hf / ^(2)+ 

+ E f <$<*r, T 4^]<F^(1)F A(T (2)>). (3) 

Here ds^ u {i) is the surface element, z = 1,2, and the double integration 
in © is performed over the minimal area surface inside the contour C. 

The Gaussian correlator (F MI/ (l)Fx CT (2)) = 1)^^(1,2) can be rewritten 
identically in terms of two scalar functions D(x) and D\(x) [3], which have 
been computed on the lattice [S| to have the exponential form D(x), D\(x) ~ 
exp(— \x\/T g ) with the gluon correlation length T g w 0.2 fm. 

As the next step one introduces the einbein variables and u; the first 
one to transform the proper times Sj, 7$ into the actual (Euclidean) times 
ti = zjj^. One has [To] 

2^) = ^, ^ ds i (D^z^) xy = constjDfi l (t t )(D 3 z^) xy . (4) 

The variable 1/ enters in the Gaussian representation of the Nambu-Goto 
form for SVmn and its stationary value vq has the physical meaning of the 
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energy density along the string. In case of several strings, as in the baryon 
case or the hybrid case, each piece of string has its own parameter i/W. 

To get rid of the path integration in (P) one can go over to the effective 
Hamiltonian using the identity 

G Q q(x,y) = (x\exp(-HT)\y) (5) 

where T is the evolution parameter corresponding to the hypersurface chosen 
for the Hamiltonian: it is the hyperplane = const in the cm. case |16j . 

The final form of the cm. Hamiltonian (apart from the spin and pertur- 
bative terms to be discussed later) for the qq case is [TBI EE] 

Ho = £M + tf , m\ , L 2 /r 2 



2u t 2 ; 2M1 - C) 2 + fi 2 ( 2 + Jo W - C)M/*)] 



+ ^r^ + r^dp. (6) 



u{(3) 

Here ( = (fix + J f3vdf3) / '(fix + fi 2 + f (3vd/3) and fit and v((3) are to be 
found from the stationary point of the Hamiltonian 

9H °\ -o 9H °\ -o (7) 

Note that H contains as input only m 1 ,m 2 and a, where are current 
masses defined at the scale 1 GeV. The further analysis is simplified by the 
observation that for L = one finds u^°> = err from (J7J) and fli = y/m 2 + p 2 , 
hence Hq becomes the usual Relativistic Quark Model (RQM) Hamiltonian 



H (L = 0) = Y, vX 2 + P 2 + ^. (8) 

i=i 

But H is not the whole story, one should take into account 3 additional 
terms H se if, H spin and H Cou i, namely, spin terms in (J3J) which produce two 
types of contributions: self-energy correction [19J 

2 Am 2 (i) An 
H self = J2—^, Am 2 = — -r/K), 77(0)^1-0.9, (9) 

•=1 2 fli 7T 

where r)(rrii) is a known function of current mass m; and spin-dependent 
interaction between quark and antiquark H sp i n |20] which is entirely de- 
scribed by the field correlators D(x), Dx(x), including also the one-gluon 
exchange part present in Dx(x). 
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Finally one should take into account gluon exchange contributions [8], 
which can be divided into the Coulomb part Hc ou i = — f^r-^, and H rad 
including space-like gluon exchanges and perturbative self-energy corrections 
(we shall systematically omit these corrections since they are small for light 
quarks to be discussed below). In addition there are gluon contributions 
which are nondiagonal in number of gluons n g and quarks (till now only 
the sector n g = was considered) and therefore mixing meson states with 
hybrids and glueballs [21]; we call these terms H m i X and refer the reader to 
21J and the cited there references for more discussion. Assembling all terms 
together one has the following total Hamiltonian in the limit of large N c and 
small T g : 

H = H + H sei f + H spin + H Co ui + H rad + H mix . (10) 

We start with Hq = Hr + H str ing where H str i ng is the term proportional 
to L 2 in (6) and Hr - all the rest terms in (6). The eigenvalues M$ of Hr 
can be given with 1% accuracy by [22] 

M 2 ^8aL + 47ra(n + ^) (11) 

where n is the radial quantum number, n = 0,1,2,... Remarkably M ~ 
4/x , and for L = n = one has /x (0,0) = 0.35 GeV for a = 0.18 GeV 2 , 
and ho is fast increasing with growing n and L. This fact partly explains 
that spin interactions become unimportant beyond L = 0, 1, 2 since they are 
proportional to dridT 2 ~ ^—^ dtidt^ (see (jHJ) and [UEH])- Thus constituent 
mass (which is actually "constituent energy") /x is "running". The validity 
of Hq as a socially accepted " constituent mass" is confirmed by its numerical 
value given above, the spin splittings of light and heavy mesons [2B] and 
by baryon magnetic moments expressed directly through Ho-, and being in 
agreement with experimental values [23]- 

We now come to the gluon- containing systems, hybrids and glueballs. 
Referring the reader to the original papers [HJ-JI] one can recapitulate the 
main results for the spectrum. In both cases the total Hamiltonian has the 
same form as in (|10jh however the contribution of corrections differs. 

For glueballs it was argued in [TT] that H © has the same form, but 
with rrii = and er — * a a dj = fcr while H sd f = due to gauge invariance. 

We now coming to the next topic of this talk: hybrids and their role in 
hadron dynamics. We start with the hybrid Hamiltonian and spectrum. This 
topic in the framework of FCM was considered in [3 EH] The Hamiltonian 
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H for hybrid looks like PQ, 0,(101 



Ht yb) = ^1 + + K + K + n* + p|±p^ + , f i _ ! + F (12) 

2 M i 2/x 2 2 2/i ^' 3 1 V ; 

Here p^, are Jacobi momenta of the 3-body system, H se if is the same 
as for meson, while H sp i n and Hcoui have different structure JU]; H str is the 
string term similar to H str ing in (10). 

The main feature of the present approach based on the Backgroung Per- 
turbation Theory (BPTh), is that valence gluon in the hybrid is situated at 
some point on the string connecting quark and antiquark, and the gluon cre- 
ates a kink on the string so that two pieces of the string move independently 
(however connected at the point of gluon). This differs strongly from the 
flux-tube model where hybrid is associated with the string excitation as a 
whole. 

Results for light and heavy exotic 1 h hybrids also given in pQ and are 
in agreement with lattice calculations. Typically an additional gluon in the 
exotic (L = 1) state "weights" 1.2^-1.5 GeV for light to heavy quarks, while 
nonexotic gluon (L = 0) brings about 1 GeV to the mass of the total qqg 
system. 



3 Hamiltonian and Fock states 

As was mentioned above the QCD Hamiltonian is introduced in correspon- 
dence with the chosen hypersurface, which defines internal coordinates 
lying inside the hypersurface, and the evolution parameter, perpendicular to 
it. Two extreme choices are frequently used, 1) the cm. coordinate system 
with the hypersurface £4 = const., which implies that all hadron constituents 
have the same (Euclidean) time coordinates X4 = const, i = l,...n, 2) the 
light-cone coordinate system, where the role of 24 and x^ is played by the 
x + ,x+ components, x + = 2fi ^p ■ 

To describe the structure of the Hamiltonian in general terms we first as- 
sume that the bound valence states exist for mesons, glueballs and baryons 
consisting of minimal number of constituents. To form the Fock tower of 
states starting with the given valence state, one can add gluons and qq pairs 
keeping the J PC assignment intact. At this point we make the basic sim- 
plifying approximation assuming that the number of colors N c is tending to 
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infinity, so that one can do for any physical quantity an expansion in powers 
of 1/N C . Recent lattice data confirm a good convergence of this expansion 
for N c = 3, 4, 6 and all quantities considered [23] (glueball mass, critical 
temperature, topological susceptibility etc.). 

Then the construction of the Fock tower is greatly simplified since any 
additional qq pair enters with the coefficient 1/N C and any additional white 
(e.g. glueball) component brings in the coefficient 1/N%. In view of this in 
the leading order of 1/N C the Fock tower is formed by only creating additional 
gluons in the system, i.e. by the hybrid excitation of the original (valence) 
system. Thus all Fock tower consists of the valence component and its hybrid 
equivalents and each line of this tower is characterized by the number n of 
added gluons. Then, the internal coordinates {£}„ describe coordinates and 
polarizations of n gluons in addition to those of valence constituents. 

We turn now to the Hamiltonian H, assuming it to be either the total 
QCD Hamiltonian Hqcd, or the effective Hamiltonian H^ e ^\ obtained from 
Hqcd by integrating out short-range degrees of freedom. We shall denote the 
diagonal elements of H, describing the dynamics of the n-th hybrid excitation 
of s-th valence state (s = rn{ff}, gg, 3g, ^{/l/a/s}) for mesons, 2-gluon and 
3-gluon glueballs and baryons respectively with /; denoting flavour of quarks) 
as H^2- For nondiagonal elements we need only the lowest order operators 

fs) (s) 

Hn n+i an d H^-in describing creation or annihilation of one additional gluon, 
viz. 

H- qqg = g J g(x, 0)a(x, 0))g(x, 0)d 3 x (13) 

H g2g = 9 -f abc J (<9 M < - d u al)ay v d*x. (14) 

As it is clear from (|13|) . (J14|) . the first operator refers to the gluon creation 
from the quark line, while the second refers to the creation of 2 gluons from 
the gluon line. In what follows we shall be mostly interested in the first 
operator, which yields dominant contribution at large energies, and physically 
describes addition of one last cross-piece to the ladder of gluon exchanges 
between quark lines, while (|T4^) corresponds in the same ladder to the a s 
renormalization graphs. 

The effective Hamiltonian in the one-hadron sector can be written as 
follows 

H = H (0) + V (15) 
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where is the diagonal matrix of operators, 

H® = {H&\H&\H&\...} (16) 

while V is the sum of operators (|13|) and (j!4|). creating and annihilating 
one gluon. In is the Hamiltonian operator for what we call the 

"n-hybrid", i.e. a bound state of the system, consisting of n gluons together 
with the particles of the valence component. In this way the n-hybrid for the 
valence p-meson is the system consisting of qq plus n gluons "sitting" on the 
string connecting q and q. 

Before applying the stationary perturbation theory in V to the Hamil- 
tonian (|15|) . one should have in mind that there are two types of excita- 
tions of the ground state valence Fock component: 1) Each of the operators 
H^,n = 0, 1, ... has infinite amount of excited states, when radial or orbital 
motion of any degree of freedom is excited, 2) in addition one can add a gluon, 
which means exciting the string and this excitations due to the operator V 
transform the n — th Fock component ip^ into V'i+i- 

The wave equation for the Fock tower ^ N {P, £} has the standard form 

H* N = (# (0) + V)* N = E N V N , (17) 

where N numerates energy eigenvalues, and £ is a set of internal quantum 
numbers in the n-hybrid, or in the integral form 

y N = - G®V* N (18) 

where is diagonal in Fock components, 

G(0, < £ >=5™b> a &w= s ->mzz (19) 

and is the eigenfunction of H^°\ 

= 40)^(0) (2Q) 

and since is diagonal, has only one Fock component, = ip n {P, {£}n) 
n — 0,1,2, and the eigenvalues E$ contain all possible excitation energies 
of the n-hybrid, with the number n of gluons in the system fixed, 



eP = E<?HP) = /P 2 + %. (21) 
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Here {k} denotes the set of quantum numbers of the excited n-hybrid. 

From (|18|) one obtains in the standard way corrections to the eigenvalues 
and eigenfunctions. 

As a first step one should specify the unperturbed functions ^f^ > intro- 
ducing the set of quantum numbers {k} defining the excited hybrid state for 
each n-hybrid Fock component ip n (P{C,} n ); we shall denote therefore: 

=ipn{k}(P,{£}n), n = 0,1,2,... (22) 

The set of functions VVffe} with all possible n and {k} is a complete set 
to be used in the expansion of the exact wave-function (Fock tower) 

®N = C l{k}^m{k}- (23) 

m{k} 

Using the orthonormality condition 

/ ^m{k}^n{ P }dT = 5 mn 5 {k}{p} (24) 

where dT implies integration over all internal coordinates and summing over 
all indices, one obtains from (17) an equation for c m ^} and En, 

C n{p}( E N - E n{p}) = Y C m{k} V n{p},m{k} (25) 
m{k} 

where we have defined 

V n { p }, m {k} = J V^{ p }^m{fc}dr. (26) 

and is the eigenvalue of the wave function component ip n {p}- 

Consider now the Fock tower built on the valence component VM K b wnere 
v can be any integer. For = 0{0} this valence component corresponds 
to the unperturbed hadron with minimal number of valence particles. For 
higher values of the Fock component ip v i K \ corresponds to the hybrid 

with v gluons which after taking into account the interaction is "dressed 
up" and acquires all other Fock components, so that the number N in (|23|) 
contains the "bare number" z/{k} as its part iV = ... (at least for small 

perturbation V). 

One can impose on ^ n the orthonormality condition 

J * N * M dT = £ C {k} c% {k} = 5 NM . (27) 

m{k} 
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Expanding in powers of V, one has 



N(u{k}) r r . N(l) . JV(2) . / OQ N 

C 44 = S m»8{kH*} + C m{k} + C m{k} + - ( 28 ) 

= + E® + ES ) + - (29) 

It is easy to see that E$ = 0, while for one obtains from (J25j) the 
standard expression 

Ml) _ Vn{p},v{K} 

C n{p} ~ p(0) p (0) • W 
^u{k} fi n{p} 

In what follows we shall be interested in the high Fock components, v+l, {k}, 

obtained by adding I gluons to the valence component v{k}. Using ([23)1 and 
(J28j) one obtains 

^,N{v{k}) _ \ - V u+ i{k},v+l-l{k!} Vv+l-l{ki},v+l-2{k 2 } 
{fcl}...{fc;} ^v{k} ~ ^u+l{k} ^u{k} ~ ^v+l-likt} 

V v+X{h} , v{K } _ +Q( yl+2y ( 31 ) 



P (0) p (0) 
U V { K } ^u+l{ki} 



Since V is proportional to g, one obtains in the perturbation series 
in powers of ot s for c N and hence for ty N (J33J). One should note that a s {Q 2 ) 
is the background coupling constant, having the property of saturation for 
positive Q 2 [HI and the background perturbation series has no Landau 
ghost pole and is defined in all Euclidean region of Q 2 . 

The estimate of the mixing between meson and hybrid was done ear- 
lier in the framework of the potential model for the meson in [27]. In [21] 
the mixing between hybrid, meson and glueball states was calculated in the 
framework of the present formalism and we shortly summarize the results. 
One must estimate the matrix element between meson and hybrid wave 
functions taking the operator V in the form of (|13|) . where the operator of 
gluon emission at the point (x, 0) can be approximated as 



x [exp(ik • x - i/it)ej i A) c A (k) + e^c^k) exp(-zk • x + ifit)] (32) 
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Omitting for simplicity all polarization vectors and spin-coupling coeffi- 
cients which are of the order of unity, one has the matrix element 

V Mh = ^f= / ^/(r)^+(0, v)d 3 r (33) 

where <pM{ r ), M VVi( r i> r 2) are meson and hybrid wave functions respectively, 
and in (}3*3*|) it is taken into account that the gluon is emitted (absorbed) from 
the quark position. 

Using realistic Gaussian approximation for the wave functions in ()33j) one 
obtains the estimate [21] 

V Mh « g ■ 0.08 GeV. (34) 

A similar estimate is obtained in J2I] for the hybrid-glueball mixing matrix 
element, while the meson-glueball mixing is second-order in (|34|). 

Hence the fist-order hybrid admixture coefficient (|3T)|) for the meson is 

n _ VMh VMh / qc .n 

Mh ~ E<$ - E<® ~ ^Mm~h (35) 



and for the ground state low-lying mesons when AMmh ~ 1 GeV it is small, 
Cmk ~ 0.1 — 0.15, yielding a 1-2% probability. However for higher states in 
the region M M > 1.5 GeV, the mass difference AM Mh of mesons and hybrids 
with the same quantum numbers can be around 200 MeV, and the mixing 
becomes extremely important, also for meson-glueball mixing, which can be 
written as 

L-MG - 2^ 7T77 A m y 6b ) 

AM Mh AM hG 

and VMh ~ VhG- A good example is given by the three fo mesons: /o(1390), /o(1500) 
and /o(1710) studied on the lattice in 

The authors [28J arrive at the following result of careful lattice studies: 

|./o(1710) >= 0.859|# > +0.302|ss > +0.413|nn >, 

|./o(1500) >= -O.128|0 > +0.908|ss > -0.399|nn >, 

|/o(1390) >= -0.495|# > +0.290|ss > +0.819\nn > . (37) 

From (|3*7j) it is clear that a strong mixing occurs between states in the 
region 1.4-1.7 GeV, however the dominant valence component in all three 
cases is clearly visible: it is (nn) for / (1390), (ss) for / (1500) and the 
two-gluon glueball (g) for / (1710). 
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4 Sea quarks in the spectrum and strong de- 



cays 

In all discussion above the dynamics in the qq system was described by the 
Wilson loop W a (C) in ([TJ. However the formalism in is correct only 

in the large N c limit, which holds for the lower part of the spectrum with 
accuracy of the order of 10%. For higher excited states with the radius ex- 
ceeding 1.3-1.4 fm one should take into account the admixture of quark pairs. 
This can be accomplished formally replacing W a {C) in by the product 
W a (C) Ylf det(m/ + D). Using the heat kernel (Fock-Feynman-Schwinger) 
representation for the determinant one has 

Re lndet(m + /)) = -lndet(m 2 + £) 2 ) = 

cfx / — exj>(-m 2 s)(Dz) xx exp(-K)W(C x ). (38) 
Jo s 

where W{C) is the Wilson loop without spin factors, present in (2), and C x 
is the closed contour beginning and ending at the point x. In this way the 
determinant can be expanded as a series in powers of number of sea-quark 
loops and the averaging yields 

(W a {C) det(m + D)> == (W CT (C))(det(m + £>)) + ^j-W x {C, C x ) + ... (39) 
where 

W 1 (C,C X ) = (W a (C)W(C x )) - (W a (C)){W(C x )). (40) 

It was shown in that the interaction of Wilson loops given by (J40|) ef- 
fectively produces holes in the original world sheet of the string in W a {C). 
It was argued in jHUj, that for unstable states with the life-time T ~ lfm 
and average radius larger than 1.4fm the holes due to the qq pairs can be 
in metastable equilibrium which can be called "the predecay state". In this 
state the linear potential between quarks starts to saturate, and the result- 
ing meson masses calculated in j30] are strongly decreased by 200-500 MeV, 
which is in good agreement with experimental data. The resulting radial 
Regge trajectories are of surprisingly good linear form, and are given in [3*T] . 

Another phenomenon, where sea quarks play the crucial role is the strong 
decay process. As it was shown in one can distinguish three possible 
mechanisms of strong decays, where nonperturbative QCD is involved: 1) 



12 



chiral mechanism, which appears after bosonization j33] 2) string break- 
ing mechanism and 3) pair creation via the intermediate hybrid or glueball 
formation. 

Taking into account valence (perturbative) gluon field a At (x), one has 
after averaging over background gluon field and bosonization the following 
effective quark-meson Lagrangian [3*2]-|34| 



J QML 



d A x / d A y{ f ipaa( x ) \[i{d-iga) + im f ] al3 5 (A) (x - y)5 fg + 



+tMsU^ } (x,y)\ i> a p(y) - 2n f [J(x,y)}~ 1 M 2 s (x,y)} (41) 

where U = exp(i<fi(x, y)^), and (f> = t a <p a / f n is the pionic field, rif is num- 
ber of flavours, and J(x, y) is the kernel proportional to the integral of field 
correlator, f,g - flavour indices, while Mg(x,y) is the effective quark mass 
operator containing the string connecting quark to the closest antiquark po- 
sition, for details and notations see [3*2^-[3^j. 

All decay mechanisms listed above are contained in the Lagrangian ()41)1 . 
e.g. the chiral mechanism obtains by expanding U(<f)) in powers of <f) a , as for 
the details and comparison with experiment, see [SB] and the literature cited 
in there. 

The nonperturbative string breaking mechanism is described in (}4*Tj) by 
the term ip + Msip, where it is essential that the mass operator Ms(x,y) = 
Mg enters at the vertex of the quark-antiquark formation (rather than the 
string mass operator in the quark propagation). The explicit computation of 
in the string breaking term was done in j22], and yields 

AL (br) = Wj* f tf x $f x ty( x y (42) 

\/7T J 



As it is clear from (|42j) one obtains the Pq mechanism with the predicted 
coefficient, which agrees with the generic phenomenologically fitted value. 

Finally, the mechanism 3) is given by the combination of the qq genera- 
tion by the perturbative gluon, due to the term ip + aip, with the subsequent 
string world sheet enveloping the gluon and quark trajectories. In this way 
one obtains the Hybrid- Mediated Decay (HMD) and - for the OZI vio- 
lating processes - the Glueball-Mediated Decay (GMD), suggested in 
The relative role and theoretical estimates of all three mechanisms was not 
elaborated in 
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5 Conclusions 



It was shown above that the Vacuum Correlator Method (VCM) is a powerful 
tool for the investigation of all effects in the QCD spectrum. In particular 
the low-lying part of the spectrum, where the string breaking is not essential, 
is well described by the Hamiltonian containing the minimal number of 
input parameters: current quark masses, string tension and a s . Here effects 
of mixing and decay are inessential (less than 10%) and the leading large 
N c approximation is valid. However for higher masses (and for high-energy 
processes) one should take into account the Fock-tower structure of hadrons, 
where in particular hybrids contribute significantly. Superficially these Fock 
towers are similar to the light-cone Fock states considered in j^Hj, however 
the (essential) difference is that in our case all Fock components are bound 
states - mesons or hybrids in the leading large N c limit, while in [35] only 
perturbative dynamics is present. 

It was demonstrated above that also the 1/N C effects, namely glueball ad- 
mixture in mesons and decay amplitudes, are easily computed in the frame- 
work of the VCM, and moreover no new parameters are introduced. However 
only first step is done in [32] for the construction of decay and production 
amplitudes. 
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